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Abstract 

We generalize the Goldberger-Wise mechanism and study the stabihty of the Crystal 
Universe models. We show that the model can be stabilized, however for configurations of 
Crystal Universe in the absence of fine-tuning, brane crystals are not equidistant, i.e. a 
" — h" pair is far away from adjacent " — h" pair, except for the fixed points of the orbifold, 
which differs from the assumptions taken in the literature. 



-•^ Mailing address in China 



Recently there has been considerable interest in studying the brane universe models. 
Randall and Sundrum (RS) [H have proposed a high dimensional picture to solve the gauge 
hierarchy problem between the Planck scale and the electroweak scale. In RS's scenario 
there are two 3-branes with opposite tension which sit at the fixed points of an jZi 
orbifold with AdS^ bulk geometry, the gravity is shown to be localized on the brane of 
positive tension (" + " brane) and the exponential warp factor in the spacetime metric 
generates a scale hierarchy on the brane of negative tension (" — " brane). If the standard 
model (SM) fields reside on the " — " brane, this leads to the resolution of the gauge hierarchy 
problem, however even though it is plausible in string theory, it is still questionable whether 
we can live on " — " brane. There are many variants of RS model proposed recently to 
avoid this problem. In Ref.0, Lykken and Randall proposed a multi-brane model, i.e. 
a " + H — " brane configuration in which the SM fields reside on intermediate " + " brane 
and a warp factor accounts for scale hierarchy. Kogan et al. in Ref.0 (see also and 

) considered " H h" multi-brane model and the Crystal Universe model (see also 

10' 111' 0' 1131) • These models provide a way of having the visible sector on a " + " brane 
with an hierarchical warp factor and interestly predict gravity be different from we expect 
not only at small scale but also at ultralarge scale. Furthermore it has been argued in [|lT| 



that this class of models give rise to many new phenomena such as neutrino mixing. Dark 
Matter which can be tested experimentally. 

In this paper we study the issue of stability of the crystal universe models. Following 
Goldberger and Wise (GW) [T^ we introduce a bulk scalar field into the models, then 
minimize the potential generated by the bulk scalar with quartic interaction localized on 
two 3-branes. We will show that the brane crystals is not equidistant in the absence of 
fine-tuning and generally a " — h" pair is far away from adjacent " — h" pair. 

To begin with, we consider a Crystal Universe Model shown in Fig.l which consists of 
n array of parallel 3-branes with " + " brane every other " — " brane in a AdS^ space with 
negative cosmological constant A. The fifth dimension y has orbifold geometry / Z2. The 
77, -f- 1 array of parallel 3-branes are located at 7/0 = 0, ?/2 ••• Un^ where ?/o = and ?/„ are 
orbifold fixed points. The action for this configuration is 



\ (1) 



S= I d^x ! dyVG {2M^R - A) - [ 

where z — 1 = 0, 1, 2, ...n, gf^V are the induced metric on the branes, Vi-i are their tensions 
and M is the 5D fundamental scale. The 5D metric ansatz that respects 4D Poincare 
invariance is given by 

ds^ = e-^''^y^r]f,^dx''dx'' - dy^, (2) 
here the warp function <j{y) is essentially a conformal factor that rescales the 4D component 



1 



of the metric. Substituting (2) into the Einstein equations we have that: 



i-1 



where k = y is effectively the bulk curvature in the region between the two adjacent 
brane crystals. There are two solution to equations (3): 



i-1 



(i) . a.(y) = (-ir+^% + ^2(-iy+i%„ (4) 

j=0 

and 

i-1 

(ii) . cj,{y) = i-iyky + J2 2(-iy%„ (5) 

j=0 

here ai{y) are warp factors between the [i — l)th brane and ith brane. For solution (i) the 
" + " brane sits on the fixed point yo = and the corresponding brane tensions are V^_i = 
(— l)*~^A/fc; for solution (ii) on the fixed point is the " — " brane and Vi_i = {—lyA/k. 

To study the stability of Crystal Universe, we introduce and couple a bulk scalar field 
to the brane crystals. This technique is a generalization of the GW mechanism, however 



the calculation invlolved in this paper will be much more complicated than that in [12 
For a bulk scalar field with mass m, 



(6) 



where Gab is the 5D metric given in (2) with cr{y) given in (4) and (5). And the boundary 
potentials of scalar field are 

S^~l = - j d'x\^i^/^{<l>'-vl,)\ (7) 

where f j_i are the vacuum expectation values of bulk scalar field in the {i — l)th brane, Ai_i 
are coupling contants. We are interested in those configurations of the bulk scalar where 
the boundary potentials are minimised. This essentially amounts to negligible dynamics 
of $ along the direction tangential to any of the 3-branes. This assumption is reasonable 
because we focus on the stability of the 3-branes system at the moment and do not study 
phenomenologcal consequence of possible coupling of the bulk scalar field $ to matter fields 
living on the branes. It, therefore, suffices to concentrate on equation of motion of $ only 
in y direction, which is 

(9J$ - 4a^(y)9^$ - = 0, (8) 
where o"-(y) = dai{y)/dy. Solution of this eqution is 

$(y) = exp {2a[{y)y)[Ai exp {a[{y)viy) + B, exp {-a[{y)uiy)]. (9) 
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In Eq.(9) Vi — •\/4 + w? /(T''^{y) is independent of y (which we denote by v in the following 
discussion), however Ai and Bi vary in the range of y^-i < y < i/i. 

To determine the coefficients Ai and B^, we require $ minimize the boundary potential. 
We firstly consider the case (i) where the " +" brane is at yo — 0, then discuss the case (ii). 

For case (i): 

A^ = C'-r'!:S^^ )Y^, (10) 



B. = ( ""'-^1^+::'^' )Yir, (11) 



where i — 2j — 1 with j — 1, 2, 3... 
and 



A = ( "-^iff'^ i-r'. (12) 



where i = 2j with j' = 1, 2, 3... 

In Eqs. (10-13), 1^ and Ri are defined as : = exp {—kyi), Ri — exp {—kyi) = Fi, and 
-Rj = exp [-k{yi - yj_i)] = Yi/Yi_i {i^l). 

Substituting A^, B^ in Eqs.(10-13) and $(y) in Eq.(9) into the action (6) and integrating 
out y give rise to a 4D effective potential V{Ri, Vi), 

A;" V(i?i, Vi) = /l(-Rl, Vo, Vi) + -^h{R2, V2, Vi) + -^[MRs: V2, Vs) 

K2 K2 

+ §/4(^4, V4, V3)] + ...+ p4^^~^ [/2j-l(-^2j-l, -^2^-2, Vsj-l) 

ri2...-n;2j_2 

^4 

+ -^f2j{R2j:V2j:V2J-l)] + ..., (14) 
^2j 

where / is defined as 

(^ + 2){R'^u - R^vf + [v-2){u- R^+'vf 
f{R, u, v) = 1-R21' ■ ^^^^ 

The effective potential in (14) is an iterated function with many variables, however it 
can be shown that a minimum exists in certain conditions. Defining that r(2j-2,2j-i) = 
V2j-2/v2j-i, r(2j,2j-i) = V2j/v2j-i, we obtaiu that 

V - V I -^2^-1/^2^+1 x2x7 

V2j-1 = V2j-1 + p4 l~ ) ^2j+l 

^2j V2j-1 
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^4 

= f2j-l{R2j-l,r(2j-2,2j-l), 1) + d4 /2i(-^2i, T(2j,2i-1), 1) 

+ %i(^)^Vi.«. (16) 

Substituting (15) into (16) V2j-i can be rewritten as 



2i-i 



(^ + 2)(ig-^._,r(2,-2,2,-l) - i?2Vl)^ + - 2)(r(2,-2,2,-l) - R'^t-if . , 



-L -n:2i-i 



(i. + 2)(i?| - i?^,r(2,,2,-i))^ + (^ - 2)(r(2,,2,-i)-i^^f)^ ~ r^2^•+l^2v1o^ 



R2j 



V2j-1 



Prom (17), one can see that for arbitrary positive values of v, Vzj-i grows as i?2j-i 1 
or as i?2j — > 1 as long as V2j-i ^ V2j-2 and V2j-i ^ V2j- These two limits correspond 
to the (2j — l)th " — " brane approaching the " + " brane ai y = 2/2j-2 and at y = y2j, 
respectively, and in these limits V2j-i is singular, i.e. V2j-i{R2j-i, R2j) ~ (1 
as R2j-i 1 and V2j-i{R2j-i, R2j) ~ (1 



Rluy > 

> as i?2j ~^ 1- This implies that the 



(2j — l)th " — " brane experiences repulsive forces exerted on it by the " + " brane of its 
either side and consequently the numbers of the branes can not be reduced. We note that 
when 1^2^-1 = ^2i-2 and/or V2j-i = V2j, the leading singularity in V2j-i is removed and the 
subleading terms in V2j-i is attractive. In this case, therefore, the less brane crystals will 
be more stable. 

Prom now on, we assume that Vi takes different numerical values in the different branes. 
To obtain the values of i?2j-i and i?2j, we minimize the effective potential Vzj-i. Por i?2j 
it satisfies the following equation: 



■■(2j,2j-l)(^2j 



where 



Since Ro 



Q 



2j 



2j+l\ 



2 (u^ + 2 -l^ + u {R%^ - 1 

,V2j+1.2 



{v + 2)R% + {v^ - 4)i?| + 2 



V 



(19) 



(20) 



i,2j ^ or 1, V2j_i — > oo which one can see from (17), this extremum represents 
the minima of V2j_i in the i?2i-dhcction. To sec whether simultaneous minima in the 
i?2i-i-direction exist, we extremize V2j_i with respect to R2i-\ and get 



(2i-2,2i-l) 



■ Jji2j-\) 



T>2-v 

2u 



2 (1 - RlU) + ^ (1 + ^U) ± ^2?'l(l - ^2--i)a/Q^] , 

(21) 
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where 

,2 



f2j{R2j,r(^2j,2j~l), 1) V2j+lV2j+l 

p4 + p4 ,,2 

2j -'^2j^2j-l 



(22) 



Note that for r(2j„2,2j-i) = ''^(2^-2 2j-i)' ^^^2j-i/f^-R2i-i < > which corresponds to a se- 
quence of saddle points. While for r(2j_2,2i-i) = ''^(2^-22^-1)' ^^^2j-i/f^-R2i-i > ^^i^ 
extremum represents the minima of V2j-i in the i?2j-i-direction. Thus in the parameter 
space where r^j2j-i) ^i^d r'^2j-2 2j-i) co-exist, the absolute minima of V2J-1 exists. There- 
fore, the minima of V{Ri,ri) exists, i.e. Crystal Universe can be stabilized. 

Having shown the possibility of stabilizing the Crystal Universe models, we discuss 
and analyze the configuration of the brane crystals when they are stabilized. Following 
GW|T^, we will also limit ourselve to the regime where e is small, e = z/ — 2^^^1. 
For a stabilized Crystal Universe model, from (21) we have 

f2j{R2j,r(2j,2j-l)A) , V2j+l f V2j+l -^2^^ fr,r,\ 
^2j ^2j ^2j-l 

Note that the two terms on the left-handed side of eq.(22) are positive, we have seperately 

f2j{R2j,r^2j,2j-l), 1) 



Rij 



V2j+1 jV2j+l 2 



< e, (24) 



< e. (25) 



R2j ^23 -\ 

In Fig. 2 we plot the allowed region of i?2j based on (23), from which one can see that for 
e = 0.01 R2j varies from 0.85 to 1, which corresponds to the (2j — l)th " — " brane very 
close to the 2jth " -|- " brane of its right side. With i?2j in this range we have r(^2j,2j-i] ~ 1, 
i.e. V2j ^ V2j-i. 

From eq. (24) we get 

^V2j+1^2^2j+l , Rjj+l ^ ^2j+1 ^2^ ^2j+3 ^2T7 (nfW 



Thus 



'V2j-1 Rij R2jR2j+2 V2j-l V2j+l 

lV2j+l N2^2j+1 



. P4 <e, (27) 

V2j-l -K2 



''2J 



p4 

^2j + l ,V2j + l.2(V2j+-i.2'r; 



-^27-^21+2 '"2j-l V2j+l' 



V2J+3 < e. (2^ 
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Given that f2j-i ~ V2j, we have from (26) 



f2j+l{R2j+l,r^2j,2j+l)A) , 
:2 ^ ^^2i, 



^(2i,2i+l) 



(29) 



1 -^2j+l 

12 d4 d4 f2j+2{R2j+2,r(2j+2,2j+l), 1) < £• (30) 

'^(2j,2j+l) -f^2j-f^2j+2 

In Fig. 3 we plot the allowed range of -R2J+1 from which we see that with i?2j in the range 
of 0.85 ~ 1, R2j+i varies from to 0.8 ~ 1 

Now we consider an additonal constraint on -R2J+1 from (27). Combining Eqs.(26), (27) 
and (29) we obtain 

R2j+l/^f2j,2j+l) ^ -^27+25 (31) 

which we plot in Fig. 4 for the allowed range of R2j+i- 

From Figs. 3 and 4 one can see that a Crystal Universe can be stabilized for a large 
parameter space of -R2j+i, however is required to be very close to 1. 

Similarly we can discuss the parameter space for solution (ii) and our results show to 
have a stabilized Crystal Universe R2j+i approachs 1, and correspondingly r^2j+i,2j) ~ 1- 
However we should point out that the absolute value of Ri can not be fixed and it depends 
mostly on the ^(i^o)- 

In summary, in this paper by explicit calculation we show that Crystal Universe can be 
stabilized by introducing a bulk scalar field to brane system. Our results differ from the 
assumptions taken in the literature for the discussion of the Crystal Universe. For exmple, 
in Refs.[B|,[0]. they have assumed that the branes are equidistant. 



We should point out that Choudhury et al. ||T3[ studied the stability of the " H h 

brane configuration and find that the " — " brane chooses to stay close to the visible " + 



brane. Taking n = 2, we recover the results of Ref. |13|. For n = 1 we agree with GW's 



results. So our results apply for general Crystal Universe models. 
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Figure 1: Crystal Universe is made up n + 1 array of " + " and " — " branes with lattice 
spacing {yi — yo), (?/2 — Vi) ■■■ [Un — Un-i) and bulk curvature k. The fifth dimension y has 
orbifold geometry /Z2, 2/0 = and |/„ are orbifold fixed points. 
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Figure 2: The allowed parameter space of R2j- The y-axis is g2j = i.e. the left-handed 

side of eq.(23); the x-axis is i?2j- The figure on the right-handed side is an amplification of 
the left in the range of i?2j — 0.85 ~ 1. 
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Figure 3: The allowed parameter space of i?2i- The y-axis is g2j = 2 — i.e. the 

''(2j,2j+l) 

left-handed side of eq.(28); the x-axis is i?2j+i- 
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Figure 4: The allowed parameter space of i?2j- The y-axis is /i2j 
left-handed side of eq.(30); the x-axis is R2j+i. 



-^^^ i.e. the 

''(2j,2j + l) 
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Figure 5: Illustration of a Crystal Universe model with " + " brane at orbifold fixed point 
1/0 = 0. The up figure ends with " + " brane and the down with " — " brane. 




Figure 6: Illustration of a Crystal Universe model with " — " brane at orbifold fixed point 
yo = 0. The up figure ends with " + " brane and the down with " — " brane. 
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